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Integrated Sachs- Wolfe effect in Cross-Correlation: 
The Observer's Manual 
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The Integrated Sachs- Wolfe (ISW) effect is a direct signature of the presence of dark energy 
in the universe, in the absence of spatial curvature. A powerful method for observing the ISW 
effect is through cross-correlation of the Cosmic Microwave Background (CMB) with a tracer of 
the matter in the low redshift universe. In this paper, we describe the dependence of the obtained 
cross-correlation signal on the geometry and other properties of a survey of the low redshift universe. 
We show that an all-sky survey with about 10 million galaxies, almost uniformly distributed within 
< z < 1 should yield a near optimal ISW detection, at ~ 5a level. In order to achieve this level 
of signal-to-noise, the systematic anisotropics in the survey must be below ~ 0.1%, on the scale of 
~ 10 degrees on the sky, while the systematic error in redshift estimates must be less than 0.05. 

Then, we argue that, while an ISW detection will not be a good way of constraining the con- 
ventional properties of dark energy, it could be a valuable means of testing alternative theories of 
gravity on large physical scales. 

PACS numbers: 98.65.Dx, 98.62. Py, 98.70.Vc, 98.80.Es 



I. INTRODUCTION 

One of the fundamental pieces of our understanding of 
modern cosmology comes from the study of anisotropics 
in the Cosmic Microwave Background (CMB). First dis- 
covered by the DMR experiment on the COBE satellite 
in the early 90's 0, the observations of CMB anisotropies 
matured through various ground-based/balloon-borne 
experiments (see for a list) until its latest climax 
by the first year data release of the observations of the 
WMAP satellite in 2003, a cosmic variance limited 
map of the CMB sky with a resolution of < 0.5 degs. 

While most of the fluctuations seen by WMAP and 
other CMB experiments were generated at z ~ 1000, low 
redshift (z < 20) physics generates additional fluctua- 
tions |3j. These secondary anisotropies can be detected 
through cross-correlating CMB measurements with large 
scale structure observations RM. Following the WMAP 
data release, various groups Q claimed a possible obser- 
vation of such correlation with various galaxy surveys, 
although at a small (2 — 3er) significance level. Most re- 
cently, claimed a correlation between WMAP maps 
and the 2MASS galaxy catalog, at both large and small 
angles. 

In this paper, we focus on such correlations on angles 
larger than a few degrees, and their only known cosmo- 
logical source, the Integrated Sachs- Wolfe (ISW) effect 
[8j . The purpose of this work is to provide comprehensive 
overview of various observational aspects and systematic 
limitations of an ISW detection in cross-correlation, in 
a flat universe. Sees. II and III will briefly review the 
physics of the expected signal and error in any ISW de- 



tection. In Sec. IV, we consider over what redshift range, 
and what angular scales, most of the ISW signal arises, 
and what the optimum signal-to-noise for an ISW de- 
tection may be. Sec. V deals with the limitations of 
realistic surveys, i.e. the Poisson noise, systematic con- 
taminations, and redshift errors. Finally, in Sec. VI, we 
discuss what we may (or may not) learn from a detec- 
tion of the ISW signal in cross-correlation, and Sec. VII 
concludes the paper. 

Throughout this paper, unless men- 
tioned otherwise, we use the flat 1st year 
WMAP+CBI+ACBAR+2dF+Ly-cv concordance cosmo- 
logical model (fi m = 0.27, fl A = 0.73, and cr 8 = 0.84 0), 
with a cosmological constant (i.e. w = —1), and a run- 
ning spectral index (n s = 0.93, and dn s /d\nk = —0.031 
at k = 0.05 h/ Mpc). However, note that our conclusions 
are not sensitive to the details of the model. 

In this paper, the 3D Fourier transform is defined as 

A k = J d 3 x A(x) exp(-ik.x), (1) 

while the cross-power spectrum of fields A and B, Pa.b, 
is defined through 

Re{A k B£,) = (27r) 3 «£(k - k')P A , B (k), (2) 

where 5^, is the 3D Dirac delta function. The auto-power 
spectrum is also defined in a similar way. 

II. THE ISW EFFECT IN 
CROSS-CORRELATION 
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The ISW effect is caused by the time variation in the 
cosmic gravitational potential, as CMB photons pass 
through it. In a flat universe, the anisotropy in CMB 
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photon temperature due to the ISW effect |8( is an integral 
over the conformal time rj 
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where $ = 9$ /9?y, and n is the unit vector along the line 
of sight. ?7ls and f?o are the conformal times at the last 
scattering and today, while the speed light is set equal to 
unity (c = 1). The linear metric is assumed to be 

ds 2 = a 2 (rf){[l + 2$(x,77)](V - [1 - 2$(x, 7y)]dx • dx}, 

(4) 

(the so-called longitudinal gauge) and 770 is the conformal 
time at the present. 

In a flat universe, in the linear regime, <I> does not 
change with time at any given comoving point for a fixed 
equation of state [3(| , and therefore observation of an 
ISW effect is an indicator of a change in the equation of 
state of the universe. Assuming that this change is due to 
an extra component in the matter content of the universe, 
the so-called dark energy, this component should have a 
negative pressure to become important at late times |lOj . 
Therefore, observation of an ISW effect in a flat universe 
is a signature of dark energy. 

The ISW effect is observed at large angular scales be- 
cause most of the power in the fluctuations of <& is at 
large scales. Additionally, the fluctuations at small an- 
gles tend to cancel out due to the integration over the 
line of sight. 

We are interested in finding the cross-correlation of 
the ISW effect with the galaxy distribution. Assuming 
Gaussian initial conditions, and full-sky coverage, differ- 
ent harmonic multipoles are statistically independent in 
the linear regime. Therefore, as the ISW effect is only 
important on large scales which are still linear, the statis- 
tical analysis is significantly simplified in harmonic space. 
For a galaxy survey with the average comoving density 
distribution n c (r) as a function of the comoving distance 
r, the Limber equation can be used to approximately 
describe the exp ected cross-correlation with the galaxy 
distribution |3lj (see Q for detailed derivation) 
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where S 2B im and Tg m are the projected survey galaxy 
overdensity and the CMB temperature in the spherical 
harmonic space. If we assume that the galaxies follow the 
matter density with constant bias, b g , (i.e. S g — b g S m ), 
then P<s>i t g(k), the 3D cross-power spectrum of $' and 
galaxy overdensity for the comoving wave-number k, can 
be related directly to P$', m . We then use the Goo Ein- 
stein equation [Tj| 

(k 2 + m 2 )>$> + m<S>' + ATTGa 2 ( Pm 5 m + PDE 5 DE ) = 0, (6) 

to relate the matter auto-power spectrum P m , m (fc) = 
P(k) to P$' t3 . Here, H = dlna/dr/ is the conformal Hub- 



ble constant, and pde and See are the average density 
and overdensity of the dark energy, respectively. 

Note that, for a cosmological constant (a ACDM cos- 
mology), Sde = 0, and Eq. JfJJ) reduces to the Poisson 
equation for k^>TL. In this case, Eq. (JSJ reduces to 
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where g is the linear growth factor of the gravitational 
potential, <f>, and g' is its derivative with respect to the 
conformal time, while z denotes the cosmological red- 
shift, which is related to the comoving distance r via 

^- = Hdr. (8) 
l + z 

For any alternative theory of dark energy, an indepen- 
dent equation for the evolution of See should be solved 
simultaneously. 



III. THE ERROR IN ISW DETECTION 

It is easy to see that the expected dispersion (see Q for 
details [22|) in the cross-correlation signal for harmonic 
multipole C' g T{i) is given by 
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where / s k y is the fraction of sky covered in the sur- 
vey, and we assumed a small cross-correlation signal, i.e. 
C 2 gT {£) < Cgg{i)C TT {t), which is the case for the ISW 
effect (the expected ISW effect is much smaller than the 
primary anisotropies, but see 01)- 

Ctt(^) is the observed CMB temperature auto-power, 
which includes both the intrinsic CMB fluctuations and 
the detector noise. As will become clear later on, since 
the ISW effect is observed at small £, the WMAP ob- 
served auto-power spectrum|l3j|. which has negligible de- 
tector noise at low £, should give the optimum power 
spectrum to use in Eq. © . We again use the Limber ap- 
proximation to obtain the projected galaxy auto-power 
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Fig. shows an example of the expected cross-power 
signal and error for a survey with 10 million galaxies with 
b g = 1, uniformly distributed between < z < 1. 



IV. THE PERFECT GALAXY SURVEY 

To obtain the optimum signal-to-noise ratio for an ISW 
detection in cross-correlation, we assume that we have an 
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FIG. 1: The expected cross-power spectrum of the ISW ef- 
fect, for an all-sky survey with b 2 g dN/dz — 10 7 , and z max = 1 



(at least approximate) redshift estimate for each galaxy 
in the survey. Then we can divide the survey into almost 
independent shells of thickness Sr, where 

r <C Sr <C r, 

where ro < 5/i _1 Mpc is the galaxy auto-correlation 
length. Using Eq. (J5J, the expected ISW cross-power 
signal for a narrow redshift bin (z, z + Sz), at multipole 
£, is 
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while Eq.@ gives the error in the cross-power spectrum 
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and Ea. HlOl) can be used to find C gg (£, z) 
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Dividing Eq. (fTT|l by Eq. lfT2]). we find the expected 
signal-to-noise ratio for the cross-correlation signal in 
multipole £, due to this shell 



6(S/N) 2 = 



C 2 gT (£) [r- 2 Sr-(2£ + l)]xiPl tm (k) 
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where k = ^ + ^ 2 . Within the approximation of indepen- 
dent shells and multipoles, (S/N) 2 is cumulative, and we 




80 



60 



40 



20 



FIG. 2: Figures show the expected (S/N) 2 distribution for 
a redshift or resolution limited full-sky ISW cross-correlation 
signal (i'max refers to the scale at which the (white) detector 
noise is equal to the true signal) . In either figures, the enclosed 
area for the region covered by a survey, multiplied by its sky 
coverage, gives the optimum (S/N) 2 for the cross-correlation 
signal. The spikiness of the distribution in ^-space is due to 
the use of the actual observed WMAP power spectrum in Eq. 
114H . Note that for partial sky coverage, the low-£ multipoles 
that are not covered by the survey should also be excluded 
from the area under the curves. 
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could simply add (or integrate over) the contribution due 
to different multipoles and shells that are included in the 
galaxy survey, and multiply it by the sky coverage, /sky, 
to obtain the optimum (S/N) 2 (in the absence of system- 
atics) for the whole survey. Fig. J5J) shows the (S/N) 2 
density distribution for hypothetical all-sky surveys with 
limited (CMB) resolution, or redshift depth, in a ACDM 
cosmology, while we assumed that the Poisson noise is 
negligible (n c — > oc). 

We see that the ISW cross-correlation signal is widely 
distributed over a redshift range between and 1.5 
(which is the era in which dark energy becomes cosmo- 
logically important), and peaks at z ~ 0.4, below which 
the detection is limited by the available volume. Almost 
all the signal is due to multipoles with I < 100 (9 > 2°), 
which implies that the WMAP 1st year all-sky tempera- 
ture map 0, with S/N > 1 on scales larger than a de- 
gree (£ < 200), captures almost all the ISW signal in the 
CMB sky, while the upcoming 2nd year WMAP data, or 
higher resolution CMB observations are unlikely to make 
a significant difference. The angular scale of the cross- 
correlation signal decreases with the depth of the survey. 
This is due to the fact that the angular correlation length 
of the galaxy distribution is smaller for a deeper (more 
distant) sample. 

For our assumed ACDM cosmological model, the to- 
tal S/N which could be achieved by a perfect survey is 
~ 7.5. Let us compare this result with previous such es- 
timates. The first estimate of S/N for a perfect survey 
is, in fact, completely consistent with 7.5a. Other 
such estimates [j| focus on specific surveys, but seem to 
be roughly consistent with the values expected from Fig. 

for those surveys. This is despite our use of Limber 
(Eq. [SJ and independent redshift bins (Eq. I19|) approx- 
imations, and demonstrates that these approximations 
are adequate for our purpose. Our approach allows us to 
visualize where in redshift-angular space one should look 
for the ISW signal. However, given a sample of galaxies in 
existence, we should emphasize that neither we, nor pre- 
vious works, provide an optimum procedure to estimate 
this signal. For a realistic galaxy survey, such proce- 
dure would depend on possible survey contaminations or 
systematics, detailed geometry, and redshift information 
available for the survey (See the next section). 

V. SUB-OPTIMAL GALAXY SURVEYS 

In this section we look into how different features of 
realistic surveys may reduce the optimum ISW signal-to- 
noise obtained in the last section. 



A. Poisson Limited Surveys 

For a realistic survey, an additional source of noise is 
Poisson fluctuations in the galaxy number density. The 
Poisson noise (the second term in the brackets in Eq. I10|) 
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FIG. 3: The distribution of (S/N) 2 for different galaxy red- 
shift distributions. For partial sky coverage, the result should 
be multiplied by /sky- 
is inversely proportional to the average number density 
of observed galaxies in the survey, which dominates the 
uncertainty in cross-correlation for a small galaxy sample. 

Figs. 111 show how the total cross-correlation 

(S/N) 2 and its distribution depend on the number of 
galaxies in an all-sky survey. Although a fixed number 
of galaxies per unit redshift is assumed, different curves 
in Fig. (PJ can be combined to obtain the signal-to-noise 
for an arbitrary redshift/bias distribution. 

Fig. Q shows that an ambitious all-sky survey with 
about 10 million galaxies (or one million clusters with 
b g ~ 3), which uniformly cover the redshift range be- 
tween to 1, can only yield a S/N of ~ 5. Although 
the Sloan Digital Sky Survey (SDSS; U\) is unlikely to 
get a S/N better than 4<r due to its incomplete redshift 
and sky coverage, future galaxy surveys like LSST and 
Pan-STARRS could achieve a cosmic variance limited de- 
tection of ISW in cross-correlation , should they cover 
almost all of sky, and have sufficiently small systematic 
errors (see Sec. IVBII . 



B. Observational Contamination 

Another feature of a realistic galaxy survey is the pres- 
ence of artificial structure in the sky maps. Such struc- 
tures may be due to observational systematics, such as 
day to day changes of weather, seeing or photometric cal- 
ibrations. Another possibility is the effective change of 
magnitude limit due to Galactic extinction, or Galactic 
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FIG. 4: The integrated (total) (S/N) 2 for a constant dN/dz 
up to Zmax, corresponding to the area under the curves in Fig. 

m 

contamination due to stars misidentified as galaxies. De- 
pending on if these contaminations are correlated with 
possible contaminations of the CMB maps, they can add 
extra random or systematic errors to the optimum error 
that may be achieved by a survey. 



FIG. 5: The graph shows the maximum allowed amplitude 
of anisotropies, Ag V3 introduced by possible systematics of 
the galaxy survey. The solid line shows the level of system- 
atics that reduces the (S/N) 2 for an ISW detection by 50%. 
The dotted line shows the level of correlated contamination 
that introduces a systematic error in an ISW detection com- 
parable to the cosmic variance. Here, we assume that the 
correlated amplitude of CMB contamination is ~ lfiK, con- 
sistent with the expected leftover Galactic contamination in 
WMAP foreground-cleaned W-band map pjjl . 



1. Uncorrelated Contaminations: Random Errors 

If the contaminations of the galaxy survey is not cor- 
related with the contaminations in the CMB map, they 
should not affect the expectation value of the cross- 
correlation function of the galaxy survey with the CMB. 
We do not expect the observational systematics of the 
galaxy/CMB surveys to be correlated. Therefore, the 
only effect will be to increase the random error through 
increasing the auto-powers of each survey (Eq. EJ . Since 
the Galactic/systematic contributions to the power spec- 
trum of the CMB are believed to be small and under 
control [T3L ITt| . we focus on the contamination of the 
galaxy survey. 

We again divide our galaxy sample into indepen- 
dent redshift bins. The harmonic components of the 
anisotropies in each redshift bin can be divided into real 
and systematic parts: 

S g (£, to; z) = S g {£, to; z) + 5f s (i, to). (15) 

While the systematic contaminations, S g ys (£,m), is also, 
in general, a function of redshift, for simplicity, we ignore 
this dependence. Of course, in practice, the level of sys- 
tematics in the redshift range that contains the bulk of 



the ISW signal is the relevant quantity (see Figs. I2l.'>ll . 
Now, the contaminated auto-power spectrum will be 

C gg (e; Zl ,z 2 ) = C gg {l; Zl )5 ZuZ2 +Clf{i), (16) 

where C gg (£,z) was defined in Eq. Q13f) . and we ignore 
the cross-term, as we do not expect any correlation be- 
tween the survey systematics and the intrinsic fluctua- 
tions in the galaxy density. Note that, unlike the pre- 
vious sections, different redshift bins are now correlated, 
and we will need to invert the full covariance matrix to 
estimate the expected signal to noise. Eq. © can be gen- 
eralized to get the covariance matrix of the cross-power 
spectrum for f s k y = 1 

C e , ZuZ2 = (AC gT (i,z 1 )AC gT (£,z 2 )) 
_ C gg (l;z u z 2 )C TT (i) 

For a low level of systematics, we can use a Taylor ex- 
pansion to invert the covariance matrix 

c 1 ~ c- 1 - c^c^cr 1 . (18) 

Now combining Eqs. I|11I12I and I16ll"5|l yields an ap- 
proximate expression for the suppression of (S/N) 2 due 
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to the presence of contaminations: 

- t (2f+ /; 9S w^(^i)c 9 T(^ 2 ) 
■J2(m +i)— 



5(S/N) 2 



C TT (£)C gg (£,z 1 )C gg (£,z 2 ) 

e^mV- (19) 



Let us define the amplitude of the systematic contam- 



inations, A g ys as 



A sys = Ml+l)C 9 T(fl 



(20) 



Again, for simplicity, we assume that A* 1 "' is indepen- 
dent of £ Of course, the relevant amplitude for A s g vs 
corresponds to the Grange that contains the bulk of the 
ISW signal, i.e. t ~ 20 (see Fig. The solid line in Fig. 
iJSJ shows the amplitude of uncorrelated contamination, 
required to reduce the (S/N) 2 by 50%. Thus, we see that 
since we need z max ^ 1 to capture most of the ISW sig- 
nal (see Sec. II V|) . the level of systematic anisotropy in 
the galaxy survey should be less than 0.1% on scales of 
£~20. 

As an example, let us assume the systematic anisotropy 
in an optical survey is due to uncertainties in the mag- 
nitude calibration. Furthermore, let us assume that we 
know the (approximate) redshifts, and in a given redshift 
bin, the luminosity function is given by the Schechter 
function [20| : 

<p(L) dL = (f>* ■ (L/L*) a exp(-L/L*)(dL/L*) (21) 

where <fi*, L*, and a are characteristic galaxy number 
density, luminosity, and slope of the luminosity function, 
respectively. Assuming that the survey is almost com- 
plete in the given redshift bin, i.e. L m i n <C L* , the sys- 
tematic anisotropy is given by: 



= (j)(L)AL min 
9 f0{L)dL 



T(l + a) 



(22) 



For a ~ —0.9 [l7|. we can ignore 1 + a in the exponent, 
and use AL/L = 0.4 ln(10)Am, to arrive at 



A sv*„ °' 4 lnl0A 



r(i + a) 



~0.lAm<10" 3 , (23) 



where Am is the uncertainty in the magnitude calibra- 
tion. Thus, we see that Am needs to be < 0.01, for the 
error to be dominated by cosmic variance. The accu- 
racy of this calibration is currently ~ 0.02 for SDSS [lH |. 
and will remain a challenge for future surveys such as 
Pan-STARRS and LSST. Note that the magnitude un- 
certainty of individual galaxies, as long as they are not 
correlated on the scale of a few degrees on the sky, can 
be significantly larger. 

Repeating this estimate for an incomplete survey, i.e. 
L m in > L* , we see that a significantly higher precision of 



magnitude measurement will be required to capture the 
ISW signal: 



Am < 



0.01 

1(r3 i5; 



for L m i n <C L* , 
for L min > L*. 



(24) 



2. Galactic Contaminations: Systematic Errors 

As the Galactic sources of contamination in the CMB 
and galaxy surveys are correlated, they may cause a sys- 
tematic over/underestimate of the cross-correlation sig- 
nal. In order to model this effect, we assume that this 
correlated signal follows the Schlegel, Finkbeiner, and 
Davis (SFD) Galactic extinction map [l!| ■ Since the bulk 
of the ISW signal comes from I ~ 20, we normalize the 
amplitude of the contaminating signals at £ = 20. The 
level of Galactic contamination in the foreground cleaned 
W-band of the WMAP satellite is expected to be 
~ 1 fiK outside the Galactic plane, and thus we take 



_ ^ _j£(£ + l)C s g f{l) 



2n 



for £ 



1 fiK, 
20. 



(25) 



while we assume C s g v g s cx cx Csfd, for the in- 

dependence of the correlated contamination. Here, Csfd 
is the angular auto-power spectrum of the SFD extinc- 
tion map. 

The dotted line in Fig. J5J shows the amplitude of the 
correlated systematic anisotropy, normalized at £ = 20, 
which yields a systematic error in the ISW signal equal to 
the cosmic variance error. Since this systematic error is 
proportional to A^ s x A g ys , the upper bound, A^ AX , 
shown in this figure decreases for larger Galactic con- 
taminations in the CMB, inversely proportional to A^ s . 
Thus, we see that since the solid curve is below the dotted 
curve, as long as Ay s < few [iK , the level of uncorre- 
lated contamination puts a more stringent constraint on 
the systematics of the galaxy survey. As we discussed in 
Sec. IV B II this requires the systematic anisotropy of the 
survey to be less than 0.1 %. 



C. Redshift errors 

Since the ISW kernel has such a broad distribution in 
redshift space, redshift errors are not expected to be a big 
limiting factor for ISW cross-correlation measurements. 

We can simulate random errors in the estimated red- 
shifts, via re- writing Eq. (|19fl for finite redshift bins of 
thickness 2Az. It turns out that even using bins as thick 
as A 2 ~ 0.5 will not decrease the signal-to-noise by more 
than 5%. 

Systematic errors in redshifts, or in redshift distribu- 
tion, may also systematically bias the estimates of ISW 
amplitude. Taking the unbiased ISW amplitude to be 
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FIG. 6: (Left) The w - Q m constraints (68% and 95% con- 
tours) based on an optimum ISW detection. (Right) The 
significance of ruling out c 3 = 0, for a quintessence model 
with given w, and c a — 1, based on an optimum ISW detec- 
tion. In both graphs, other cosmological parameters are kept 
constant. 



one, the ISW amplitude, biased by the systematic red- 
shift error Az is given by 



1 + 
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ISW 
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C gT (£,z + Az)/C gT {£,z) 
AC 2 gT (£,z)/C 2 gT (£,z) 
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which, for small Az, yields 
| (SIS W | 



(26) 



ISW 

„ Az-dC*(£,z)/dz 

E ' 



Ac 2 aT (e,z) 
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E 



^ 3 Az, 



(27) 



where we used Eqs. (I11I13|) to estimate the numerical 
value of the above expression. Since the optimum signal- 
to-noise for an ISW detection is ~ 7, the redshift system- 
atic errors may not dominate the errors if 



A;. 



< 



1 



3x7 



0.05. 



(28) 



VI. 



WHAT DOES ISW TELL US ABOUT 
COSMOLOGY? 



Let us study the optimum constraints that an ISW 
cross-correlation detection can give us about cosmology. 
For a nominal concordance cosmology Co, the expected 
significance level for ruling out the cosmology C, Xj is 
given by 



x 2 = 



E 



[C gT (e,z;C)-C gT (e,z;Co)]* 
AC 2 gT (l,z;C ) 



(29) 



where C g T and AC g T are defined in Eqs. i|llll3[l . Note 
that the bias factor, b gi is cancelled from the numerator 



and the denominator, and the optimum significance level 
only depends on the fluctuations in matter density and 
gravitational potential. 

Fig. |JBJ shows the optimum constraints that an ISW 
detection may yield on some of the properties of dark 
energy, i.e. its density, equation of state, and speed of 
sound ]22J . While such constraints are already compara- 
ble to the current bounds on these parameters pA l2lll23l . 
future observations of CMB and large scale structure |2J| 
will significantly improve these bounds. However, ISW 
constraints could still be used to test possible systematics 
that may affect the accuracy of future measurements. 

It has been claimed |2^ that the cross-correlation ob- 
servations could be used to break the degeneracies in the 
so-called doomsday cosmological scenarios, which involve 
a varying equation of state. However, given the cosmic 
variance limited measurements of the ISW effect (Eq. 
I29fl . it does not seem that such constraints can signifi- 
cantly improve the expected combined constraints from 
the upcoming Planck and SNAP satellite experiments 

mm 

A more intriguing application of an ISW detection is 
the possibility of testing our theories of matter/gravity 
on large scales 28] . While the consistency of the power 
spectrum of 2dF and SDSS galaxy surveys [2!| with 
the WMAP power spectrum confirms the consistency 
of the ACDM concordance cosmology at scales of k > 
0.1 h Mpc -1 , Fig. |0| shows that an ISW detection can do 
the same at scales of k ~ 0.003 — 0.03 h Mpc -1 . There- 
fore, current 0, and future observations of an ISW 
effect may confirm the consistency of our cosmological 
theories at the largest physical scales that they have ever 
been tested. While the cross-correlation statistics are al- 
most free of systematic bias, the auto-correlation is often 
dominated by survey systematics on such scales. 



VII. CONCLUSIONS 

In this paper, we study different aspects of the corre- 
lation between a galaxy survey and CMB sky, due to the 
ISW effect in a flat universe. The main source of noise 
is contamination by the primary CMB anisotropies. We 
see that, given this noise, most of the signal comes from 
£ ~ 20, and z ~ 0.4 with negligible contribution from 
£ > 100 and z > 1.5. 

While the maximum significance level for an ISW de- 
tection is ~ 7.5c for a concordance ACDM cosmology, 
an all-sky survey with about 10 million galaxies within 
< z < 1 should yield an almost perfect ISW detection, 
at the ~ 5<r level. 

We show that, in order to achieve a cosmic- variance 
limited ISW detection, the systematic anisotropies of the 
galaxy survey must be below 0.1%, on the scale of ~ 10 
degrees on the sky. Within this level of systematics, 
and given the current level of Galactic contamination in 
WMAP CMB maps (~ fiK), there should be a negligible 
over/underestimate of the ISW signal due to the system- 
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significantly reduce the significance of an ISW detection, 




the systematic errors in redshift estimates need be less 
than 0.05, in order to achieve a cosmic-variance limited 
measurement. 

It turns out that, due to the large noise induced by 
the primary anisotropies, the optimum constraints on the 
properties of dark energy from an ISW detection, are al- 
ready comparable to the current accuracies, and will be 
outdone by future observations. However, the simplic- 
ity of the linear physics involved in the ISW effect, and 
the fact that the cross-correlation statistics are not bi- 
ased by the systematics of CMB or galaxy surveys, makes 
ISW detection a useful indicator of possible systematics 
in more accurate methods. 

Finally, we point out that the detection of the ISW 
effect provides a unique test of our concordance cosmo- 
logical model on the largest physical scales that it has 
ever been tested. 
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This statement is strictly valid only for wavelengths much 

larger than the Hubble radius, or if the speed of sound 

vanishes. 

The error in the Limber approximation, in the form ex- 
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pressed here, is about 10% at £ — 2 and drops as £~ 2 . 
[32] If the survey does not cover the whole sky, the nearby 
multipoles in the harmonic space will not be independent. 
However, if we bin a large number of multipoles together, 



we can use the increased standard deviation of Eq. © for 
each multipole and treat them as independent to estimate 
the error of each multipole bin. 



